ASYMPTOTIC BEHAVIOR FOR A NONLOCAL 
CONVECTION-DIFFUSION EQUATION 



LIVIU I. IGNAT, TATIANA 1. IGNAT, DENISA STANCU-DUMITRU 

Abstract. In this paper we consider a nonlocal convection diffusion equation and obtain 
the first term in the long time behavior of the solutions. The results of this paper are 
obtained by using a scaling argument and a new compactness argument that is adapted 
to the nonlocal character of the considered problem. 



1. Introduction 

The aim of this paper is to analyze the long time behavior of the solutions of the following 

system 

f Ut = J*u-u + G* - X G M'', t > 0, 

(1-1) 

I w(0) = y^, 

where J and G are smooth positive functions with mass one. Moreover we assume that J 
is radially symmetric. 

The well-posedness of the model presented here has been analyzed in [T3l Th. 1.1]. 
For any g > 1 and G L^(]R'^) fl L°°(R°') there exists a unique global solution u G 
C([0, oo), L\W^) n L°°(M'^)) satisfying 

Since J and G have mass one the mass conservation property holds 

u{t,x)dx= / Lp(x)dx. 

Moreover as proved in [T3l Th. 1.4] the solutions decay similar to the classical heat equation: 
for any 1 < p < oo the following holds: 

(1-2) ||'u(t)||ip(Kd) < C{p,d, ||v9||ii(B;d), ||v?||L-(Md))(^ + l)"2(^~p). 

This decay property has been obtained in [13] by the so-called Fourier Splitting method 
EZ] and in a more general setting in [TB]. When p = 2 a. similar argument has been 
also used in As far as the authors know, the case p = oo in fll.2p is open. 

In the case when the nonlinear term is supercritical, i.e. q > 1 + 1/d, the first term 
in the asymptotic behavior has been analyzed in [13]. There the main idea was that the 
nonlinear part decays faster than the linear semigroup and then the first term in the long 
time behavior is given by the linear semigroup. This has been already observed in [9] in 
the case of the classical convection-diffusion system. 
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The aim of this paper is to give an answer to the critical case q = 1 + 1/ N even though 
we give a proof that both treats the critical and super-critical case. The subcritical case 
l<g<l + l/A^is still open. The method we employ is the so-called four step method 
that consists in the analysis of some rescaled orbits {u\{t)}. We refer to [32] for a review 
of the method in the case of the porous medium equation. 

We consider two important quantities 



A 



J{z)\z\^dz and B = {Bi, . . . , Bd), Bj 



G{z)zjdz,j 



d. 



The main result of this paper is the following one. 

Theorem 1.1. Let 1 < p < oo. For any Lp E L^{W'') fl L 
(11. ip satisfies 



^) the solution u of system 



:i.3) 



limt2(i p) 



uit) - t-'"'f\ 



LP 



where the profile fm is the smooth solution of the equation 



-Af^ - ■ Vfm = pm - aB ■ V(|/„ 



with J^afr, 



m where m is the mass of the initial data ip and 

0, g>l 



a 



1 

d- 



Next, we say a few words about the above asymptotic profile 

U{t,x) 



When q>l + l/doiB = Oi^^ the asymptotic profile is given by the heat kernel. When 
q = 1 + 2 ^^'i B ^ Oi^rf, U is the unique solution of the following equation 

Ut = AAU - B ■V{\U\'>-'^U), a;GM^t>0, 
f/(0) = m5o. 



:i.4) 



The well-posedness of this system has been analyzed in [10] in the one- dimensional case 
and in [11] the multi-dimensional case. It has been proved in [1] that the profile fm is of 
constant sign and decays exponentially to zero as |x| — oo. 

We remark that in the case of the symmetric function G, i.e. G{z) = G{—z), the solution 
of (II. ip converges to the heat kernel since in this case B vanishes. When B ^ we obtain 
in the limit the solutions of the viscous convection-diffusion equation. Along the paper 
we will consider the case of nonnegative initial data, so nonnegative solutions of system 
(ll.ip . The case of sign-change solutions could also be analyzed with small modifications of 
the proof (see [10] for a rigorous treatment of the critical case for the convection-diffusion 
equation). 
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In the linear case, i.e. Ut = J * u — u, the asymptotic behavior has been obtained in 
[Hj by means of Fourier analysis techniques and in [T7] by scaling methods. In [T7] the 
scaling argument works since it is applied to the smooth part of the solution. Refined 
asymptotics have been obtained in [m [16] . We also recall here [301 El] where a scaling 
method is used for equations of the type Ut = J * u — u — . There the authors obtain 
barriers for W and its derivatives, W being the smooth part of the solution of the linear 
equation Ut = J * u — u. Once these barriers are obtained the authors split the solution 
of the nonlinear problem in a way that permits to obtain uniform Holder estimates and 
then compactness. The method developed here is more flexible in the sense that it uses 
only energy estimates that involve the linear part of the equation and the good sign of the 
nonlinearity. 

In the local case, i.e. Ut = Au + a ■ V {\u\'^~^u) , the same analysis has been performed 
in a series of papers. In [5] the case q > 1 + 1/d is treated and the results in the critical 
case have been obtained by a careful space-time change of variables and using weighted 
Sobolev spaces. The sub-critical case is more difficult and the one- dimensional case has 
been considered in [10] . The extension to higher dimensions has been obtained in [11] and 

In contrast with the analysis in [9] here we assume that the initial data belong to the space 
L\R'^)r]L°°(E.'^). This assumption is necessary since even in the linear case Ut = J *u — u 
a lack of smoothing effect is present. More precisely the solutions of the linear model are 
as regular as the initial data. In the case of the heat equation with initial data in L^(M'^) 
the solution at any positive time belongs to any L^(M'^) space and this type of gain of 
integrability can also be proved for the nonlinear convection-diffusion [9]. 

We recall some similar models to those analyzed here. In [21] the author considers a 
one-dimensional model that is nonlocal in the diffusive part Ut = J * u — u + uux with 
J = e~'^' and he proves that its solutions converge to the ones of Bourger's equation with 
Dirac delta initial data. However the key tool used there, an Oleinik estimate d^uit) < 1/t 
in V{M.) is not available in our model. The methods used here can be adapted to analyze 
similar models but with nonlinearities of the type {u'^)x, > 2, [23]. In these cases, entropy 
conditions in the sense of Kruzkov [2D] should be imposed on weak solutions in order to 
have a well-posed problem. This does not appear in our model since the nonlinearity does 
not involve derivatives. 

The models considered here could be related with the ones considered in [8] where a 
scalar conservation law with a diffusion- type source of the type ut + V ■ f{u) = APgU is 
analyzed. There Pg is essentially given by Pgui^) ~ (1 + \^\'^)~'^u{C) ^-iid even more general 
models are considered. However, in order to obtain the long time behavior of the solutions, 
the authors assume that the initial data belong to some spaces where is large 

enough. The analysis of these models by our methods remains to be considered in future 
papers. 
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Similar nonlocal models have been introduced recently in [7] where the nonlocal convec- 
tive term takes the form 



where is an odd function. The possible application of the methods introduced here 
remains to be analyzed. The main difficulty in these models is even from the beginning the 
global existence of the solutions. Some models when the convection is nonlocal has been 
considered previously in [TB], Uf = u^x+G^u'^—u'^ , g > 2 and [IH] , Ut = Uxx+{u''~^{K*u))x, 
q = 2. The main difficulty in obtaining the asymptotic behavior for similar models where 
the convection is dominant, i.e. 1 < g < 2, is to obtain an entropy estimate. If the entropy 
inequality can be avoided in the critical case it seems to be crucial for the uniqueness of 
the solutions of the limit equation in the sub-critical case. However, we refer to where 
the asymptotic behavior of systems of the type Ut = Au — dy{\u\'^~'^u) with q subcritical 
is obtained without entropy estimates but rather with a kinetic formulation that allows 
to use some compactness arguments previously employed in the case of multidimensional 
scalar conservation laws [22]. The possible application of these kinetic methods to the case 
on nonlocal diffusion and/or convection remains to be analyzed in the future. 

The paper is organized as follows. In Section |2] we review a few compactness arguments 
known to be useful in the analysis of time evolution problems. We obtain a nonlocal 
compactness argument that guarantees the convergence in spaces of the form Lf^^((0, oo) x 
Mf^). Our results extend the ones obtained in [UlMlll]. The argument is similar to the three 
space method given by Aubin-Lions Lemma [29]. Once the compactness tool is obtained, 
in Section [3] we prove Theorem II. the main result of this paper. 



In this section we review a few classical compactness tools and give some results that 
will allow us to prove the main result of this paper. 

Now we recall some results given in [22] about the characterization of compact sets in 
L^(0, T, B) where B is a. Banach space and 1 < p < oo. 

Theorem 2.1 ([2S], Th. 1). Let 7 C Lp{0,T,B). 7 is relatively compact m 1^(0, T,B) 
for 1 < p < oo, or C(0, T, B) for p = oo if and only if 

(1) < / f{t)dt, f E 7? is relatively compact in B for all < ti < t2 < T . 



(2) \\Thf - /|Up(o,r-h,B) -^0 as h^O uniformly for / G J". 

The following criterion is also given. 

Theorem 2.2 ([21], Th. 5). Let us consider three Banach spaces X ^ B --^ Y where 
X ^ B is compact. Assume 1 < p < oo and 

i) T is hounded in L^(0,T, X), 

\Thf - /||LP(o,T-h,Y) as /i uniformly for / G J". 




2. Compactness Tools 
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Then T is relatively compact in Lp{0, T, B) (and in C(0, T,B) if p = oo ). 

The last criterion is obtained mainly by using Theorem 12.11 and the following inequality 
that follows from the fact that X is compactly embedded in B: for any 5 > there exists 
77(e) > such that 

(2.5) \\u\\B<e\\u\\x + vii^)\\u\\Y, \/ueX. 

In the nonlocal setting we will obtain a similar inequality in Lemma 12.21 

Now we recall some compactness results that have been proved in the nonlocal context 
[1], [3] and in a more general setting in [21] . 

Theorem 2.3 ([3], Th. 6.11, p. 126). Let I < p < 00 and n C R be an open set. Let 
be a nonnegative smooth radial function with compact support, non identically 
zero, and pn{x) = n'^p{nx). Let {/n}n>i be a bounded sequence in L^iVL) such that 



(2.6) nF / / p„(x - y)\fn{x) - fn{y)\''dxdy < M. 

Jn Jn 

The following hold: 

1. If {/n}n>i is weakly convergent in L^iVt) to f then f G W^'P{Q) for p > 1 and 
f G BV{Q) forp = 1. Moreover 

\\Vf\\L.in)<C{Q,p)M. 

2. Assuming that Q is a smooth bounded domain in M.'^ and p{x) > p{y) if \x\ < \y\ then 
{fn}n>i is relatively compact in L^{Q). 

We point out that the assumption on the compact support of function p could be re- 
moved. In fact once we have estimate (12.61) for p we also have this estimate for any other 
compactly supported function p with p < p. 

The above results hold under more general assumptions on the weights {pn}n>i and on 
a bounded domain Q in M"' with Lipshitz boundary. As proved in |2H Th. 1.2] we can 
assume that {pn}n>i is a sequence of radially symmetric functions in L^(M'') satisfying 



Pn > 0, a.e. in 



(2.7) 



and that 




j^'iPnix) = 1, V n > 1, 
lim^^oo 



^"^"^ ^^■\fn{x)-fr.{y)\Pdxdy<M. 



'nJn F — 1/1 

Then the results in Theorem 12.31 remain true in dimension d > 2. In dimension d = 1 some 
technical assumptions have to be assumed [211 Th. 1.3]. Choosing pn{x) = n{n\x\yp{nx) 
with p radial and decreasing, these technical assumptions hold and we obtain the results 
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in the second part of Theorem I2.3[ We also recall that under the above conditions on p„ 
a Poincare inequality holds [211 Th. 1.1] 



fn I f"- 



<C(p,fi,{pJ) / / Pf'' l\ Ux) - UvWdxdy. 
L^i^) JnJn \x-y\^ 



In view of this inequality the boundedness of {/n}n>i in L^{Q) is guaranteed by (12. 6p if we 
assume that {/n}n>i is bounded in L^[Q) and Q has finite measure. 

The main compactness tool that will be used in this paper is the following one. 

Theorem 2.4. Let 1 < p < oo and Q CM.'^ be an open set. Let p : M'^ — M &e a nonnegative 
smooth radial function with compact support, non identically zero, and Pn{x) = n'^p{nx). 
Let {fn}n>i be a sequence of functions in L^((0,T) x Q) such that 

(2.8) r f i/.r < M 

Jo Jn 

and 

(2.9) [ [ [ Pn{x-y)\fn{t,x)- fn{t,y)\^dxdydt<M. 



n Jn 



1- If{fn}n>i is weakly convergent in Lp{{0,T) x to f then f e LP{{0,T),W^'P{n)) 
forp>l andfe L\{0,T), BV{Q)) forp=l. 

2. Let p > 1. Assuming that fl is a smooth bounded domain in M'^, p(x) > p{y) if 
\x\ < \y\ and that 

(2.10) \\dtfn\\Lp{io,T),w-i,pm < M 

then {fn}n>i is relatively compact in L^{{0,T) x Q). 

Remark 1. Extensions to mixed type space norms of the type Lp((0, T), L^(n)) could also 
be obtained by adapting the estimates in this paper. The possibility of obtaining more 
general nonlocal compactness tools as in Theorem \2.2\ remains to be analyzed. In fl2.10p 
for technical reasons we considered the space W~^'''^{Q) but we believe that the results still 
hold by replacing W~^''^{Q) with any space Y such that L^{Q) ^ Y continuously. 

Proof. Using the same arguments an in the proof of Theorem 12.31 (see [3], Ch. 6, p. 128]) 
we obtain the results in the first part. 

We now prove the second part of the theorem by following the ideas in [29] but taking into 
account the particular estimate (12.91) . From now on, in order to simplify the presentation, 
we assume that p is a smooth radially symmetric function supported in the unit ball, 
non-identically zero and satisfying p{x) > p{y) if |x| > \y\. 

Step. I. Compactness in Lp((0, T), iy~^'^(fi)). We now check the hypotheses in 
Theorem 12.11 Let us choose < ti < t2 < T and set 

/•t2 

9n{x) = / fn{s,x)ds. 
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Estimate (12.91) gives us that 

n"" [ [ Pnix - y)\gnix) - gnivWdxdy < MT^-\ 
Jn Jn 

Theorem 12.31 applied to sequence {gn}n>i shows that there exists g G W^'^{Q) such that, 
up to a subsequence, (?„—;■ (7 in L^(f2) so in W-^'P{n). Estimate (^AO\i shows that the 
second requirement in Theorem 12.11 is also satisfied. Hence {fn}n>i is relatively compact 
in LP{{0,T),W-^'P{n)). 

Step. II. Compactness in Lp((0, T), L^^^{fi)). Since {/n}n>i is bounded in Lp{{0, T) x 
VL) then up to a subsequence {/n}n>i weakly converges to some function / in -L^((0, T) x 
The first part of Theorem [Ml guarantees that / G Lp((0,T), W^^p{SI)). 

We now use the strong convergence in -L^((0, T), H^~^'^(f2)) obtained in Step I, estimate 
(12. 9p and the fact that / G L^((0,T), ^^^'^(fi)) to prove that up to a subsequence, {/„}«>! 
strongly convergences to / in Lp((0,T) x VL). 

To simplify the presentation we will always denote the subsequence by {/«}• Also, when 
possible, we will not write all the constants in inequalities of the type / < Cg using instead 
f ^9- 

In the following we prove that for any C f2 such that d{VL' ^ dVt) > 0, {/„}n>i is 
relatively compact in L^((0,T) x Q'). From now on for a set O we will denote 

Op = O + supp(p) = {x + 9, X E Q,6 E supp p}. 

The following two Lemmas will be very useful in our analysis. Their proof will be given 
later. 

Lemma 2.1. Let Q be and open set o/M'^. For any 1 < p < 00 there exists a positive 
constant C{p,Q,p) such that the following inequality 

(2.11) nP [ [ Pn{x-y)\u{x)-u{y)\Pdxdy<C{p,n,p) I \Vu\p 

Jn Jn Jn 

holds for all n > and u G ^^^'^(fi). 

Lemma 2.2. Let Q be a bounded domain and x ^ Cl{Q). There exists a positive constant 
C = C{Q, Xi P^P) such that for every e G (0, 1) the following inequality 

(2.12) C [ \xu\P<enP [ [ p^(x - y)\u{x) - uiy)\Pdxdy + e [ \u\p + -\\ur^_, 

Jn Jnp„ Jnp„ Jn,^ ^ 

holds for all ue^Ip > 1 and for all u G Lp(R"'). 

Remark 2. In the right hand side of inequality (12.121) we have and the W~^'P{yL)- 
norm. We believe that some improvement in (I2.12p can be done by allowing the norm of 
the last term to be in some space Y with Lp{Q) Y and replacing correspondingly. 
The extension of Lemma \2.2\ to general spaces Y will enlarge the class of nonlocal problems 
where the scaling arguments used in this paper can be applied. 
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Let US &x Q' G Q and choose a smooth function x compactly supported in Q such that 
X = 1 in fi'. We choose Nq large enough such that fl'^^^ C Q for all n > Nq. 

Applying Lemma 12.21 with g = fn — f and to the set Q' we have for any e > and 
n > e~^/P that 

(2-13) ||x^llip(f7') ^ ^'^^ / / Pn{x-y)\g{x)-g{y)\Pdxdy + e I |^|^ + 

Jn' Jn' Jn' ^ 



Pn Pn pTL 



We integrate the above inequality on the time interval [0,T] and obtain that 

f [ xVn - mxdt <e7f r [ [ p„(x - y)\{fn - f){x) - ifn ' fMl'dxdy 
Jo Jn' Jo Jn', Jn'„ 



Pn Pn 

T 



e 



<FnP 



T 

Pnix - y)\fnix) - fn{y)\^dxdy 



+ enP [ [ [ pn{x-y)\f{x)- f{y)\Pdxdy 
Jo Jn'^„Jn'^„ 

^{\\fn\\LP{{0,T)xn'pJ + ll/lll,P((0,r)xn;:,,j) ~^ 5 l''^" ~ / II Vl/-i.P(n)^^- 

Since for n > ma,xNo,e~^^P we have that il'p^ C we use estimates f l2.8p . fl2.9p . Lemma 
Oand the fact that / G Lp((0,T), W^^P{n))\o obtain that 

£ j^^ xVn - ffdxdt < e(2MT + ||/rL.((o,T),VKi..(n))) + ^ /^^ " f\\w-^'nn)dt- 
Using Step I, up to a subsequence, we obtain that for any 5 G (0, 1) 

limsup [ I \fn- fl^dxdt <e{2MT ^ 

n.-5>oo Jo Jn' 



P \ 
LP{{0,T),W^'P{n)))- 



Then fn strongly converges to / in L^((0,T) x Q'). Applying a standard diagonalisation 
procedure we can extract a subsequence, denoted again by {fn}n>i, such that /„—;■/ in 

Step. III. Compactness in ^^((0, T), We now use the following result in pH 

Lemma 5.1, Lemma 7.2]. For a positive number r > we set 

^Ir := {x E : d{x, OQ) > r}. 

Lemma 2.3. Let Q be a bounded Lipschitz domain of R"'. There exist constants tq > 
depending on Q and on p and Ci, C2 (depending on p, Q and d) so that the following holds: 
given < r < Vq we can find Nq > 1 such that 



(2.14) /|^?r<Ci/ \g\P + C2rPnP [ f pn{x - y)\g{x) - g{y)\Pdxdy 

Jn Jnr Jn Jn 



NONLOCAL CONVECTION-DIFFUSION 9 
for every g G Lp{Q) and n > Nq. 

We apply the above Lemma with g = fn — f and integrate the resulted inequality on 
the time interval (0,T). Thus 

^ [ \fn-ff< r [ \fn-ff + rVr[ 

Jn Jo Jur JO 

< r [ \fn-fr + rV r [ [ prXx - y)\Ux) - fn{y)?dxdy 
Jo JUr JO Jn Jn 

+ r^n^ r [ [ pr,{x-y)\f{x)-f{y)\^dxdy. 
JO Jn Jn 

Using estimate fl2.9p and Lemma 12.11 we get 

r / l/n - /r < r / l/n - /r + r^M + r [ l V/^. 
Jo Jn Jo Jnr Jo Jn 

Since fn^fiTL LP{{0, T), Lf^^(f2)) we can let n — )■ oo and then for any r G (0, Tq) we have 

hmsup r I \f^-f\Pdx<T^(M+ f I |V/r 
n->-oo Jo Jn ^ Jo Jn 

This implies that, up to a subsequence, fn ^ f L^{{0,T), Lp{Q)) and the proof of 
Theorem 12.41 is now finished. □ 

Proof of Lemma \2.1\ We first consider the case when f2 = W^. By scaling, it is sufficient 
to consider the case A = 1. Since 

(2.15) u{x) — u{y) = I {x — y) ■ \/u{y + s{x — y))ds 



we get that 

p{x — y)\u{x) — u{y)\^dxdy < fj p{x — y)\x — y\^ \\/u{y + s{x — y))\'^dsdxdy 

Jo 

p{z)\z\P [ \Vu\P. 

JK^ 

In the case of a bounded domain fl we first extend u to M'^ such that || Vu||2,p(iHd^ < 
C(f2)|| Vu||LP(n). Then we have 



'^^ Pnix - y)\u{x) - u{y)\^dxdy < pn{x - y)\u{x) - u{y)\^dxdy 

J Jnxn J JR2d 

<c{p,p) [ <cip,n,p) [ \vuf. 

Jw Jn 
The proof of Lemma 12.11 is now complete. □ 

The rest of this subsection is devoted to the proof of Lemma 12.21 In order to give its 
proof we need some auxiliary Lemmas. 
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Lemma 2.4. Let 1 < p < oo. There exists a positive constant C = C{p,p, d) such that for 
every e G (0, 1) the following inequality 



(2.16) C\\u\ 



< e 



Pn{x-y)\u{x)-u{jj)\'''dxdy+\\u\\\^^^^i^ +s ^\\ 



holds for all ne^/P > 1 and for all u G L^iW^). 

Before starting the proof of this Lemma a few comments are needed. The case p = 2 is 
reduced after applying the Fourier transform to the following inequality 



(2.17) 



C{p)<en'[l-p{^)) + l 



Using that p is a smooth radially symmetric function we obtain that its Fourier transform 
decays at infinity and moreover, 1 — p{^) ^ \^\^ for ^ — 0. This shows the existence of two 
positive constants ci and C2 such that 



(2.18) 



This property implies that inequality f l2.17p holds for all n > 5^^/^. 
The local version of inequality fl2.16p is the following one 

(2.19) 



< 



S\\U\\\Y1,; 



= e\\{I - A)i/2«|U.(i,.) + e-'\\{I - A)-i/2«||^,( 

We remark that when p ^ 2 this inequality is not a consequence of a duality argument since 
the dual of W^'P{R'^) is W'^'^'iR'^). Inequality fl219ll holds by proving that, depending on 
the Fourier localization of u, its L^-norm is controlled by one of the two terms in the right 
hand side of f l2.19p . In fact, using classical multiplier arguments (see [121 Ch. 5]) 

\\u\\LPiR^) < eW - A)^/^n||ip(ud), suppn C : 1^1 > e"^/^} 

and 

\\u\\lpiwi) < eW - ^y^^^uWLPm^), suppu C : 1^1 < e'^^^}. 



Proof of Lemma 2.4 Let us first make a change of variable to avoid the presence of Pn{x) 
n'^pinx). Estimate (12.161) is equivalent to the following one 



(2.20) 



Lv{ 



< e 



n' 




p{x — y)\u{x) — u{y)Ydxdy + 

+e-i(/-n2A)-i/2^r 



LP 



We use a decomposition of u that has already been used in [TB]. Let us choose rj G C^(M'^) 
with 

/ rj = I and |V?7| + |?7| < p. 

This choice of r] can be always done if p is positive in some open set. We write 

u = V + w, V = rj * u, w = u — V. 
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\v\\LP{Rd) < C(?7)||u||ip(Kd), ||u'||LP(R<i) < C(?7)||u||ip(]Rd) 



We now emphasize some important properties of v and w. First of all observe that both 
of them have the L^-norm controlled by the L^-norm of u: 

(2.21) 

and moreover 

||^||lp(r<') < ||^'||lp(R'*) + 
Since the mass of rj is one we have the following representation for w: 



w{x) = 

Holder's inequality gives us that 



r]{x - y){u{x) - u{y))dy. 



(2.22) 



IP < 



p/p' 



\ri{x — y) I \u{x) — u{y) \^dxdy 



<C{ri,p) / p{x -y)\u{x) -u{y)\Pdxdy. 
In the case of v, since J^^ dxjV = 0, j = 1, . . . ,d, we write its gradient as 

Vf = V?7 * M = / V7]{x — y){u{x) — u{y))dy. 
Thus the same argument as before gives us that 

(2.23) / \Vv\P<C{r]) [ [ \Vri{x - y)\\u{x) - u{y)\Pdxdy 

<C{r],p) / p{x -y)\u{x) -u{y)\Pdxdy. 
We now prove estimate (I2.20p . In view of (12.221) for en^ > 1 we have that 

(2.24) / \w\P<enP I p{x - y)\u{x) - u{y)\Pdxdy. 



We claim that v satisfies the following inequality for all e G (0, 1) and en^ > 1 

(2.25) ||^|U.(M.) < e\\{I - n^/^f/\\\Ln^^) + ^"1^ - n^^Y^'M. 

Estimates (jMll), (l2:23|) and (|2:25|) imply that 



< 



< 



+ e-'\\iI-n^A)-'/^v\ 



LP 



< 



\u\P + nP 
\u\P + nP 



J{x - y)\u{x) - u{yYdxdy\ + e~^\\{I - n'A)-'/\r] * n)|U.(K.) 



J{x - y)\u{x) - u{y)\Pdxdy + e~^\\{I - n^A)-^^\\\ 



I LP 



Taking into account the above estimate and estimate f l2.24p for w, we obtain that f l2.20p 
holds. It remains to prove that fl2.25p holds. Writing explicitly the terms in the right hand 
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side of fl2.25p we reduce it to the case n = 1. In this case inequahty fl2.25p is exactly 
estimate fl2.19p in the local setting. □ 

Lemma 2.5. Let Q be a smooth bounded domain o/M'^ and p G (l,oo). For any smooth 
function x supported in Q there exists a positive constant C (x) such that 

(2.26) llxwllw-i.p(Rd) < C{x)\\u\\w-^,p{n)- 

Proof. We consider the case of the smooth function u. The general case follows by density. 
By the definition of the space W~^'P{M.'^) there exists a sequence G W^''^ (M"') with 
llv'nllvKi.p'(Rd) ^ 1 such that 



Since x has the support included in f2, we have x^n ^ W^'^ (il) and 
W^'^^WwyiU) - \\x\\w-'-°°in)\\Vn\\w^y (Rd) < C{x)- 

Hence 

J ^XWn = j ux^n < \\u\\w-^•p{n)\\x^n\\y^,l,p' < C{x)\\u\\w-^-pin)- 
Letting n — )■ oo we obtain the desired estimate. □ 



Lemma 2.6. Let p be a radial function with compact support, p(0) ^ 0, Q be a domain 
in M'^ and 1 < p < oo. For any x G there exists a positive constant C = C{x,p, ^) 

such that the following inequality 



(2.27) CnP / / p„(x - y)\{xu){x) - {xu){y)\'dxdy 

<n^ j j p„(x - y)\u{x) - u{y)\^'dxdy + ( j p{z)\z\A 
holds for any n > and any u G L'p{W^). 

Proof. Let us first observe that since p is radially symmetric and p(0) 7^ we have 



supp pn = — supp p. 
n 



For x ^ Vtp^ and ?/ G we have that |?/ — x| > 1/n and then p„(x — y) = 0. li y ^Vt then 
x{x) = x{y) = 0. Similar things hold if we interchange x and y. Hence 



(2.28) nP I I pn{x~y)\{xu){x)-{xu){y)\Pdxdy 

Pn{x - y)\{xu){x) - {xu){y)\^dxdy. 



Using the following identity 

{xu){x) - {xu){y) = x{x){u{x) - u{y)) + u{y){x{x) - x{y)) 
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we obtain that 

(2.29) f I Prix - v)\[xn)[x) - {xu){y)\'dxdy 



'^T^^WxWl'^m I j pn{x - y)\u{x) - u{y)\Pdxdy 

^Pn ^ ^Pn 



+ nP / pnix-y)\u{y)\P\xix)-x{y)\^dxdy. 
Using identity ( I2.15P for x it follows that 
(2.30) [ [ pn{x - y)\u{y)\P\x{x) - x{y)\''dxdy 



I / Pn{.x - y)\u{y)\^\x - y\^ I \{Vx)iy + s{x - y))fdsdxdy 
Jnp„ Jnp„ Jo 

< llxIki^-fRd) / p{z)\z\Pdz / \u{y)\Pdy. 



Putting together estimates (12:281) , fl2:29|) and (^M) we infer the desired estimate (jMlD • □ 

Proof of Lemma \2.^ From Lemma 12.41 we know that for any e G (0, 1) and ne^/P > 1 the 
following inequality holds for all v G Lp{M.'^) : 

IMlpiR'i) ^ ^1^^ [ [ Pn(a;-y)K'(a;)-i;(?/)|Prfxd?/ + ^||i;||^_i,p,jj,. 
We now localize the above inequality by applying it to v = xu- Thus 

\xu\^<£nP I I Pn{x-y)\{xu){x)~{xu){yWdxdy + -\\xu\\l^r-l.p(Kd^ + £\\xu\\lp(Rd)■ 
jRd jRd e ^ > ' 

By Lemma 12.51 and Lemma 12.61 we deduce that 

Xu\P<enPI I pn{x - y)\u{x) - u{y)\Pdxdy + e I \u\^ + -\\u\\w--^.p{n) 
and the proof is finished. □ 

3. Proof of the main result 

Before starting the proof of Theorem 11.11 we need some preliminary results that will be 
used along the proof. 

3.1. Preliminaries. In the following we denote 

Jx{x) = A^J(Ax), Gxix) = A'^G(Ax), G{z) = G{-z), Gx{x) = A'^G(Ax). 
In our nonlocal context the key compactness result is given by the following proposition. 
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Proposition 3.1. Let {fn}n>i be a sequence in {{0 , T) , L"^ [M.'^)) such that 

(3-31) |l/n.||L2((o,T)xR'*) < M, 

(3.32) ^^111 Jn{x-y){fn{x)-fn{y)fdxdy<M 

Jo JR'i Jr'' 

and 

(3-33) \\dtfn\\L'2{{0,T),H-'>-iRd)) < M. 

Then there exists a function f G L^{{0,T), H^{M.'^)) such that, up to a subsequence, 
(3.34) zn LL((0,T)xM'^). 

Proof. We apply the first step in Tlieorem 12.41 to sequence {/n}n>i and to ^2 = M"'. As- 
sumptions ( I3.3ip and fl3.32p guarantee tlie existence of a function / G L^{{0,T), H^(M.'^)) 
such that /„ weakly converges to / in L^((0,T) x M'^). The strong convergence in f l3.34p 
follows from the second step of Theorem 12.41 □ 

The following Lemmas will be used along the proof of Theorem 11.11 

Lemma 3.1. The following integration by parts identities hold 
(3.35) 

/ (J* $ - $)(x)^f(x)rfx = / $(x)(J * ^' - ^')(x)rfa; 

1 

~ ~2 

and 



J{x - y){^{x) - $(?/))(^(x) - ^!{y))dxdy 



(3.36) / {G*<^ -<^){x)^!{x)dx = / ^{x){G*^! -^){x)dx. 

Proof. Use Fubini's theorem and in the first case the fact that J{—z) = J{z). □ 

Lemma 3.2. For any p G [l,C)o] there exist two positive constants C{p,J) and C{p,G) 
such that 

(3.37) ||A2(J,*^-^)||i,(K.) <C(p, J)p2^|U.(K.) 
and 

(3.38) II A(Ga * ^ - V-) II LP(R^) < C(p, G) II V^|Up(m.) 
hold for allX>0 and G G'^iR'^). 

Proof. We treat the cases p = 1 and p = oo since the other cases follow by interpolation. 
Taylor expansion up to the second order gives us that for any x,y eW^ the following holds 

ipiy) — 'ip{x) = Wip{x){y — x) + / {1 — s){y — x)D^il}{x + s{y — x)){y — x^ds. 
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After a change of variables we have 

A2(J,*^-V;)(x) = A'^+2 f j^xix-y))[ijiy)-^ix)]dy = \' [ Jiz)Uix - ^) - ijix))dz 



= X^ f J{z) [ - f ■ V^{x) + 4 /" (1 - s)zD^ij{x - ^)z'ds 
jRd I A A Jo A . 

Since J is radially symmetric we have 

(3.39) / J{z)zj = for all j = 
and 

(3.40) / J{z)zjZk = for all 1 < j ^ A; < d. 
Those identities give us that 

(3.41) * ^ - ^){x) = V Al -s) [ J{z)zjZkj^-^{x - '-^)ds 

j^iio V dxjdxk A 

and then for p G {1, oo} inequality f l3.37p holds with C(J) = | J^a J{z)\z\'^dz. 
In the case of the second estimate (13.381) we use the identity: 

'ijj{y) — 'ijj{x) = / {y — x) ■ 'V'iIj{x + s{y — x))ds 



dz. 







(3.42) 



and apply the same arguments as in the first case. □ 

3.2. Proof of Theorem II. IL We consider the family {u\{t)}\yo defined by 

ux{t,x) = X'^u{XH,Xx). 
It follows that Ma is a solution of the following rescaled equation 

Mt = X\Jx * ux - ux) + A'^(^-'?)+2(^^ *ul-ul), X e t > 0, 

ux{0,x) = ipxix), 

where (px{x) = A°'v^(Ax). 

The proof of Theorem 11.11 is divided into four steps. 

Step I. Estimates on the rescaled solutions ux- We recall [13], Theorem 1.4] that 
solution u of system (11.11) satisfies for any p E [1, oo) and t > the following estimate 

(3.43) ||'u(t)||ip(Md) < C{p, \\ip\\Li(M.d), ||v9||ioo(iRd))(t + 1)^2(1-1). 

In the sequel we will denote by C a constant that may change from line to line, may 
depend on ||y9||j;^i(]Kd) and ||iy9||^oo(]Kd) but it is independent of the scaling parameter A. In 
the following lemmas the constant M will depend on ||<y5||^i(]{j£i) and ||(p||£^oo(iRd). We will not 
make explicit this dependence unless this is necessary. 
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Lemma 3.3. For any < ti < t2 < oo and p G [l,oo) there exists a positive constant 
M = M{ti,p) such that 

\\ux\\L^((ti,t2),LP{R'i)) < M 

holds for all X > 0. 

Proof. Using estimate fl3.43p and the fact that the rescaled solutions satisfy 

\\ux{t)\\LP(Rd) = A''*^^~^)||n(A^t)||ip(Kd), 
we deduce that for any j9 G [1, oo) and t > the following inequality holds for all A > 0: 

(3.44) \Mt)\\LnR^) < c (^p^j < 

□ 

Lemma 3.4. For any < ti < t2 < oo there exists a positive constant M = M{ti) such 
that the following inequality 

f f f Jx{x-y){ux{t,x)-ux{t,y)f dxdydt<M 

holds for all X > 0. 

Proof. Multiplying fl3.42p by ux and integrating over M'^ we get 
(3.45) 

li\\Mt)\\l2»d)= f X\Jx*ux-ux)ux{t)dx+ ! X''^^~'^'^+\Gx*ul-ul)ux{t)dx. 



'^dt 

Using that Gx has mass one the last term in the above identity is negative. Indeed, 
{Gx*ul){t,x)ux{t,x) dx = / Gx{x - y)ul{t,y)ux{t,x)dxdy 

<f f Gxix-y)(^ul''\t,y) + ^ul+\t,x))dxdy= f ul+\t,x)dx. 
Next, integrating (13.451) over the interval (^1,^2) and using identity (I3.35P we obtain 
lkA(^2)|li2md) + A^ /" f [ Jx{x -y){ux{t,x) -ux{t,y)y dx dy dt < \\ux{ti)\\l2 



Using inequality (13.440 in the case p = 2 we conclude that 

/" [ [ Jx{x -y){ux{t,x) -ux{t,y)f dx dy dt < Gt^^ 

Jti J]R<* JR'i 

and the proof finishes. □ 

Lemma 3.5. For any < ti < t2 < 00 there exists a positive constant M = M{ti) such 
that 

lh'iA,t||L2((ti,t2),H-i(Rd)) < M 

holds for all X > 1 . 
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Proof. Multiplying by e ^^(M'^), integrating over M*^ and using Lemma [3. II we get 

ux,t{t, x)i;{x) dx= [ \^{Jx *ux- ux)ij{x) dx + [ X'^^^-'^^^^G x *ul- ul) iP{x) dx 

2 

+ [ X^^^~''^-^^{Gx*^-^)ul{t,x)dx, 



J\{x - y){i'{x) - ^{y)){ux{t, x) - ux{t, y)) dxdy 



where Gx{x) = Gx{—x). Using Cauchy's inequality, the fact that A > 1 and g > 1 + l/d 
we get 

ux,t{^^^)^i^) < [— / Jxix - y){ij{x) - i){y)ydxdy] 

\^ r r o , ^ \V2 

dy 

1/2 



X (— / / Jx{x - y){ux{t,x) - uxit,y)f dxdy) 
+ \\X{Gx*^p-^p)\\L^iR'^)( f Mt,x)\"' dx 



Applying Lemma [27T] to Jx and ip, Lemma [32] to Gx and estimate (13 .44^ to ux we deduce 

that 

(3.46) 

ux,tit, x)ip{x) dx 

d 

<||^||^.(M.)[(A2 / / - y){ux{t,x) - ux{t,y)Y dxdy)'^' + t-i^'^-'^' 
Thus 

lkA,tWII?f-i(Rrf) ^ [ [ Jxix-y)iuxit,x)-uxit,y)Y dxdy + r-2^^'^"'\ 

Integrating this inequality on the time interval (^1,^2) and then applying Lemma [3.41 we 
obtain the desired result. □ 

Step II. Compactness in ^^^^^((0, 00), L-'^(M'')). We first establish the compactness in 
L^^Q^((0, 00) X Mf^). Using estimates on the tail of {ux} we will obtain strong convergence 
in LL((0,oo),LHM'^)). 

Lemma 1X51 and Lemma 1X51 give us that {ux} is uniformly bounded in L;^^((0, 00), L^^^(M.'^)) 
and {dfUx} is uniformly bounded in Lf^^{{0, 00), H~^{Q)) for any bounded domain Q of 
Mf^. Taking into account that L'^[Q) is compactly embedded in H~^{Q) for any e > 0, 
and H~^{Q) is continuously embedded in H'^{Q) for < e < 1, by classical compact- 
ness arguments ([29], Corollary 4, p. 85) we deduce that {ux} is relatively compact in 
C{[ti,t2], H~'^{Q)) for all < ti < ^2 and < e < 1. Extracting a subsequence we get 



18 L. I. IGNAT 

Using estimate (I3.44p we obtain that for each fixed t > 0, the family {ux^{t)}n>i is 
uniformly bounded in L^^^(R°'). Then any subsequence {^Afc„(^)}n>i weakly convergent 
should converge to U{t). Indeed, if u\^^(t) ^ v in L^(f2) then u\^^{t) f in V'{Q) and 
hence v = U(t). This fact shows that for every t > and p G [1, oo) we have 

ux„{t)-U{t) in LIM')- 

The uniform bound in ( ]3.44p of {ux{t)} transfers to U(t). Hence, the limit point U 
belongs to L^^((0, oo), //^(M*^)) for all 1 < p < oo and moreover we get that 

(3.47) ||t/(^)|Up(K.) <liminf ||nA(t)|Up(M^) < -j^, V t > 0. 

In the particular case p = 1, for any t > we obtain that ||f/(t)||j;^i(]jd) < m, m being the 
mass of the initial data ip. 

Let us now prove the strong convergence in Ll^^{{0, oo) xM'^). Lemma |33| Lemma [3l4l and 
Lemma [375] show that for any < ti < ^2 < oo there exists M = M(ti, \\'^\\L^(Md^, ||v^||l°°(r<*)) 
such that 

(3.48) \\ux\\Lmt^,t2)xRd) < M, 

(3.49) ^'/V / Jx{x-y){ux{t,x)-ux{t,y)fdxdydt<M 

Jti JW JR'* 

and 

(3-50) \\ux,t\\L2i{tut2),H~HRi)) <M. 

We apply Proposition 13.11 to family {ma}a>o and time interval (^1,^2)- We obtain that 
there exists a function v G L'^{{ti,t2), H^{M.'^)) such that, up to a subsequence, 

ux-^v in L\{t,,t2);LlM''))- 

The previous analysis shows that v = U. Thus U G ^^^^((0, 00), i/^M"')) nLi„^((0, 00) x W^) 
and 

ux^U in L;i„^((0,oo) X R"'). 

We now prove that in fact ux strongly converges to U in L;^^^((0, 00), L^(]R'')). Using a 
standard diagonal argument the compactness in L,i„^((0, 00), L\W^)) is reduced to the fact 
that for any < ti < ^2 < C)0 the following holds 

I't2 

(3.51) / '^Ux{t)'^L^{\x\>R)dt — 7- as i? — )■ 00, uniformly in A > 1. 

This follows from the Lemma below since the initial data belongs to L^(]R^). 

Lemma 3.6. There exists a constant C = C{J, \\'p\\L^[Rd), ||v5||L°°(Rd)) such that the follow- 
ing inequality 



t t^/^ 
^{x)dx + C{— + —) 

x\>2R ^|x|>_R -n- -n- 



(3.52) / ux{t,x)dx < I 

J\x\>2R J\i 

holds for any t > and R > 0, uniformly on X > 1. 



NONLOCAL CONVECTION-DIFFUSION 19 

Proof. Let G C°°(M'^) be such that < < 1 and satisfies iIj{x) = for |x| < 1 and 
ip{x) = 1 for |x| > 2. We put iPr{x) = ip{x/R). We multiply equation fl3.42p by ipR and 
integrating by parts we obtain 

ux{t,x)ijjR{x)dx - ifx{x)il)R{x)dx=\^ I / ux{s,x){Jx* ^Pr - i^R)dxds 



Jo JR'i 

_^y(l^,)+2 I" I ul{s,x){Gx*tpR-4jR){x)dxds. 
Jo JR'i 

We now use Lemma [3.21 with p = oo, the fact that 

ll^^(^i?)llL-(iR'') = R~^\\D^'ip\\L--iR'i) and || VV^ij||ioo(Kd) = i?~^||V^||icx>(]Rd) 
and the conservation of the (]R'^)-norm of ux to find that 
(3.53) 

ux{t, x)ipR{x)dx < / (px{x)'ijjR{x)dx + C{J)R~'^\\D'^ijj\\i^oc(^d^ / / ux{s,x)dsdx 

f I ui^(^s,x)dxds 
Jo Jr'' 

< I ^x{x)dx + C{J)R'~^D^ij\\L^^^d)t\\ip\\Li^^d) 

J\x\>R 

+ A^(i-'')+i/2-iV^|Uoc(Kd) f [ ul{s,x)dsdx. 

Jo Jr'' 

To estimate the last term in the above inequality we use the decay of ux as given by (I3.44p 
and obtain that 

Jo JR'i Jo (l + A^s)^" ^0 {l + s)^~ 

Since for any q > 1 + ^ 

ds 



and 



lim X ^ , , - 

_T , ds , 2x f 0' ^ > 1 + i 

lim X I -T, — TV = lim -r, — — = < 

(1 + .)^ + \<oo, g = l + i, 



we find that 

yd(l-q) 



^ [ [ ul{s,x)dxds<Ct^'^. 
Jo Jr'^ 

Going back to (13.531) . using that A > 1 and tp{x) = 1 for |a;| > 2 we get 

r t r I ^1/2 

ux{t,x)dx< ^ix)dx + Ci— + —)< ipix)dx + C{— + —) 

x\>2R J\x\>\R -n- J\x\>R ^ ^ 
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and the proof of the Lemma is finished. □ 

Lemma [3T6] shows that u\ ^ U in Ll^^{{0, oo), L^(M'')). This result also shows that for 
a.e. t > we have 

(3.54) \\ux{t) - f/(t)||Li(iRd) as X ^ oo. 

This fact will be used in Step IV to obtain the main convergence result of this paper. 

Step III. Passing to the limit. Using the results obtained in the previous step we 
can pass to the weak limit in the equation involving ux- Let us choose < r < t. For any 
test function ip e C^(M'^) we multiply equation fl3.42p by ip and we integrate on (r, t) x M"'. 
We get 



u\{t,x)ip{x)dx — / ux{t, x)i/j{x)dx 

X\Jx*ux-ux)ij{x)dxds + X'^'^^^''^+^ [ [ {Gx*ul-ul)ij{x)dxds 

X\Jx* - 4j)ux{s,x)dxds + X'^'-^-'^^^'^ [ [ {Gx* - ^p)ul{s,x)dxds. 
First of all observe that since for any t > 0, Ux{t) U{t) in L^^^{M.'^) we have 

ux{t, x)^jj{x)dx — / ux{t, x)^jj{x)dx ^ / U{t,x)il){x)dx — I U{t, x)^lj{x)dx. 

JR'i JR'i JRd- 

Using identity fl3.4ip and the Lebesgue dominated convergence theorem we obtain that 



2 ^,J=1 



X\Jx *^-^){x)^-y ^ I J{z)ZiZj^;^-^ = AA^(X) 



dxidxj 



where 



A = l [ J{z)\z\^dz. 

2 jRd 

Since Ma — > f/ in L^{{T,t) x M^) we obtain by using the Lebesgue theorem that 



X'^{Jx * ^ — ip)ux{s, x)dxds A / / A'ipU{s, x)dxds. 



For the term involving G we prove that 
(3.55) 

yi[i-q)+2 I I (^Gx*i'-i')ul{s,x)dxds 



0, ?>i + i 

J^J^,B-V^ix)U''is,x)dxds, q = l + i 



where 



5 = (5i,...,5d), B,= [ G{z)z,dz. 

JRd 
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When q > 1 + 1/d we use Lemma 13.21 for G and estimate (13.441) on the L^- norms of ux to 
get 

Jt JR'i 

< A'^(i-^)+iA(GA*^-^)IUo.(K.) r / ul{s,x)dxds 

* ds 



< 



C(G)||V7/;|U^(i,.) A 



(1 + A2s)l(«"i) 
< C(G)||V^^|Uoo(K.) A^"'^(^-i) ^0 as A ^ oo. 

In the case q = 1 + 1/d we prove that 

(3.56) [ [ X{Gx*iIj -'^){x)ul{s,x)dxds [ [ B ■V'^{x)U'^{s,x)dxds. 

Jt Jr^ Jt Jr'' 

Observe that 

X{Gx*i'-i^){x) = X [ G{z)U{x-^)-i){x)]dz = X I G{z)U{x + ^) - i){x))dz 

Jw \ X / J^d \ X / 

= X [ G'(z) [^V^(x) + 4 [ (l - s)zDmx + ^)z'd.s\dz. 
jRd LA A Jo A J 

Hence 

(3.57) \X{Gx *i,-i,){x)-B- V^(x)| < X-'WD^i^h^^^dy 
For any p E [1, oo) we have that 

\\ux{s)\\LP{Rd) < G{p,t), V s > r. 

Thus, up to a subsequence, 

Since ux ^ U a.e. on (r, t) x R'^ we conclude that x = U'^. Using now (I3.57P we obtain 
that fl336|) holds. 

All the above convergences show that U satisfies 



U{t, x)tlj{x)dx — / U{t, x)ip{x)dx 

= A / U{s,x)Aip{x)dxds + a I I U''{s,x)B ■V^{x)dxds, 



where a = \\iq = l + l/d and a = for g > 1 + 1/d. Thus, when g>l + l/(iori?7^ Oi^^, 
f/ is a distributional solution of the heat equation. When q = 1 + 1/ d and B ^ Oi^^, U is 
a distributional solution of the equation: ut = AAu — B ■ V{u'^). 
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Step IV. Identification of the initial data. Let us choose r = in the previous 
step. Then for any G C^(M'^) we get 



ux{t, x)t(j{x)dx — / ux{0, x)ilj{x)dx 

< ||D^V'l|L°°{Rd) / / ux{s,x)dxds + X'^^'^~''^^'^\\Dtlj\\L^(^d^ / u{{s,x)dxds 
Jo Jr'' Jo Jr'' 



ds 



fl + s)^~ 



Letting A — oo we get 

U (t, x)ilj{x)dx — m'i/'(0) 



where m is the mass of the initial data (/?. This show that for any G C, 



2(jn,d\ 



lim / U{t, x)t{j{x)dx = rmjj{0). 

no jj^d 

We want to show that this is true for any smooth bounded function and then U(t) — )■ mSo 
in the weak sense of nonnegative measures in M.'^. 

Let us now choose ip a bounded smooth function. For any e > we choose ipe G C'^iJS.'^) 
such that \\ip — ipeWioo^^d-^ < e. Then 



U (t, x)ipix) — mip{0) 



< 



U{t,x){^lj{x) -i/je{x))dx +m|?^(0) -^e(0)| -\ 
< 2em + \\ip,\\w2,^(^^d){t + t^/^). 
Thus there exists to = ^o(^) such that for all t G (0,to) the following holds 



U(t, x)ips{x)dx — mips{0) 



U (t, x)ip{x) — rmp^O) 



< Aem. 



This shows that U (t) goes to m6o as t — in the sense of measures. 

In conclusion the limit point U satisfies U G L^^ii^, oo),L\R'^))nLl^{{0, oo), H\R'^)). 
When q>l + l/doTB = Oi^, f/ is a solution of the heat equation with m6o initial data. 
When q = 1 + 1/d and B ^ Oi^^, ?7 is a solution of the equation: 



(3.58) 



ut = AAu - B ■ xeW^,t>Q, 
n(0) = m^o- 
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Since for any r > we have U{t) G -/^^(M*^) classical results on parabolic equations show 
that for any r > 

U e C((r, oo),L\W^)) n L°°((r, oo) x R'^). 

Using the fact that the heat system as well as system fl3.58p have a unique solution (see [TO] . 
[11] for complete details) then the full sequence {ma}, not only a subsequence, converges 
to U. 

Step V. The asymptotic behavior. We recall that from Step II we have 

||nA(l) - f/(l)||Li(Kd) ^ asA^cx). 
After setting t = and using the self-similar form of U {t, x) 

we obtain 

lim ||iz(t)-f/(t)|Ui(i,.) = 0. 

This is exactly fll.31) in the case p = 1. The general case, 1 < p < oo, follows immediately 
since 

2p-2 

The proof of Theorem 11.11 is now completed. 
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